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Abstract 



We discuss the scaling behaviour of different fermion actions in dynamical 
simulations of the 2-dimensional massive Schwinger model. We have chosen 
Wilson, hypercube, twisted mass and overlap fermion actions. As physical 
observables, the pion mass and the scalar condensate are computed for the 
above mentioned actions at a number of coupling values and fermion masses. 
We also discuss possibilities to simulate overlap fermions dynamically avoiding 
problems with low-lying eigenvalues of the overlap kernel. 
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1 Introduction 



Besides the interest in the 2-dimensional Schwinger model pp in its own right as a 
quantum field theory, it can be considered as a test laboratory for new theoretical 
concepts and ideas that aim at eventual applications in more demanding situations 
such as lattice QCD. In particular, for numerical simulations the lattice Schwinger 
model is most suitable to perform test studies since the computations are much 
cheaper than in four dimensions and precise results at many parameter values can 
be obtained, see e.g. refs. [21 El HI E] for a selection of recent work in the Schwinger 
model. In this paper, we want to address the scaling properties of a number of 
fermion actions for Nf = 2 flavours of dynamical fermions. To this end, we will 
compare standard Wilson j^j, Wilson twisted mass 0|B], hypercube [5] and overlap 
fermions i n their approach to the continuum limit. We will use throughout the 
paper the Wilson plaquette gauge action with a coupling (3 = l/a 2 e 2 , with a the 
lattice spacing and e the dimensionful coupling. 

Each of the above mentioned fermion actions has certain advantages and are 
used in present simulations of lattice QCD. Understanding the scaling properties 
when using different lattice actions and to check the expected scaling behaviour as a 
function of the lattice spacing a is certainly one of the most important questions in 
lattice calculations. However, if we think of dynamical fermion simulations in lattice 
QCD a scaling analysis is, at least nowadays, far too computer time consuming to be 
addressed, see e.g. ref. [TT] and ref. [T2| for recent estimates of the simulation costs 
in lattice QCD. On the other hand, for the 2-dimensional Schwinger model, such 
simulations are perfectly possible and give important insight into the properties of 
the above mentioned actions. 

In order to study the scaling behaviour, we will fix the scaling variable z = 
(m/v^) 2 ^ 3 ) where mj is the fermion mass in lattice units. We have chosen z = 
0.2, 0.4, 0.8. The fermion mass is determined from the PCAC relation where we 
employ local as well as conserved currents. At each of these fixed values of z we 
compute the pseudo scalar mass and the scalar condensate and follow its behaviour 
with decreasing value of the lattice spacing. Performing finally a continuum limit 
of our results allows us to compare to analytical predictions that are available from 
approximations of the massive Schwinger model which cannot be solved exactly. 

The paper is organized as follows. In section 2 we give the definition of the 
Wilson, hypercube, maximally twisted mass and overlap fermion operators which 
we have employed in this work. Section 3 is devoted to the observables we have used 
and provides a description of the numerical simulations, in particular our attempts to 
simulate overlap fermions dynamically. Section 4 contains our results and section 5 
the conclusions. 
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2 Lattice fermions 



In this section we give the definitions of the different kind of fermions we have used, 
i.e. Wilson, hypercube, Wilson twisted mass and overlap fermions. Throughout the 
work we have employed only the Wilson plaquette gauge action with gauge fields 
U n ,n £ U(l), n denoting a lattice point and \i — 1, 2 a direction of the 2-dimensional 
lattice 

So\U]=PY,l( U r + U *p) W 
p 

with Up denoting the plaquette 

Up = U nA U n+l2 Ul +U <2 , (2) 

where 1 and 2 denote shifts in direction 1 or 2 respectively. The coupling multiplying 
the plaquette gauge action is /? = 1/g 2 . Denoting by e the physical coupling, 
dimensions are introduced by (3 = l/a 2 e 2 . We restrict ourselves to Nf = 2 flavours 
of dynamical fermions. 

2.1 Wilson fermions 

As a first action that can serve as a kind of benchmark action, we have chosen 
standard Wilson fermions fOJ with the Wilson operator D-w given by 

= ( m + 2r)5 n _ m - - J^[(r ~ ^n)U n ,^n,m-p, + + ^)£ / rL,;A,m+/i] • (3) 

Here the sum goes over the two directions \i — 1 and \i — 2 and we use the standard 
Pauli-matrices o" M ,yU = 1,2,3 with 03 = diag(l, — 1), n,m denote lattice points 
with space and time coordinates (x, t) and we suppress the Dirac indices. The 
Wilson parameter r is chosen to be r = 1 throughout the paper. The Wilson action 
is expected to lead to large discretization errors in physical observables linear in 
the lattice spacing a and hence approaches the continuum limit rather slowly. We 
will compare results obtained with the Wilson action to corresponding results from 
actions where these 0(a) lattice artefacts are expected to vanish and changed to an 
0(a 2 ) behaviour. 

2.2 Hypercube fermions 

Perfect actions jHj are completely free of lattice artefacts. However, since such actions 
cannot be realized in practice, truncated versions are used which are expected to 
inherit many of the properties of the perfect actions. For this work, we have chosen 
a particular ansatz, the hypercube action ^H]- I n this approach the interaction 
among fields on the lattice is extended from a purely nearest neighbour interaction 
to a hypercube of the lattice. 
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A general ansatz for the hypercube operator is 

D hyp = Yl P™-n a H + X m-n 1 , (4) 
/' 

where p^-n an d A m _ n are real parameters, the values of which have to be opti- 
mized according to some criterion, e.g. by demanding an improved continuum limit 
behaviour of physical observables. 

Furthermore, symmetry requirements of the action leads, in the case of the 
Schwinger model, to a reduction of these couplings to only 5 free parameters that 
will depend on the gauge coupling constant g and the bare fermion mass m hyp . The 
values of p and A are determined for free fermions, setting g = and then they are 
taken over to the interacting case. As an alternative, the values of p and A may 
be determined by requiring that the violations of the Ginsparg- Wilson relation are 
minimized. In this work we will, however, only work with the "scaling optimized" 
set of parameters as given in ref. ^3] since we are interested mainly in the scaling 
behaviour. We give in the following table the values of the parameters that we have 
used in this work. 



Ao 


+ 1.5 + m hyp 


Ai 


-0.25 


A 2 


-0.125 


p U) 


+0.334 


p (D 


+0.083 



The full hypercube Dirac operator is given by 
D hyp = A + # (1) + # (2) 

2 

{H W T m = ^ (u n>IJ ,8 n>m -jx{\x + cr M p (1) ) + Ul_- itfl 5 n!rn+[i (X 1 - <v> (1) )) 

(H i2) ) nm = \ ((U n ,iU n+u + U n , 2 U n+%l ) 5 nim „i„ 2 (A 2 + + a 2 )p^) 

+(U n ,iUl +u%2 + Ul_ %2 U n _ %l ) ^, m _ i+2 (A 2 + (a, - a 2 )pW) (5) 
+(Ul_ u U n _ u + U n , 2 Ul +2 _ u ) 5„, m+i _ 2 (A 2 + ( a\ + o 2 )p^) 

2.3 Twisted mass fermions 

The twisted mass formulation of lattice fermions has been introduced originally to 
regulate small, unphysical eigenvalues of the Wilson lattice Dirac operator jjj. In 
order to keep an 0(a) improvement, the twisted mass setup has been first developed 
in the 0(a) Symanzik improved theory. However, it then was realized that by a 
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careful tuning of the parameters of the Wilson twisted mass action, an automatic, 
full 0(a)-improvement can be reached, leading to lattice discretization errors that 
appear only in 0(a 2 ) and hence allow for a much accelerated continuum limit [H|. 
The Wilson twisted mass formulation has received a lot of attention recently and a 
number of tests in the quenched [T^l EH El El E] and partly also for full dynamical 
fermions [201 HH 1221 122] has been performed. 

In order to introduce twisted mass fermions, let us start with the continuum, 
Euclidean action, 



S[tp, ip} = J d xipia^Df, + m + ifi f a 5 T 3 )ip, (6) 

where T3 is the third Pauli matrix, acting in flavour space and /i/ represents the 
twisted mass parameter. The transformation 

$ = ^e 1 ^ (7) 
leaves this form of the action invariant with rotated mass parameters, 

m! = m cos to + fi f sin u 

fi'f =— m sin uj + fif cos uj , (8) 
with the rotation ("twist") angle uj, 

tancj = ^. (9) 
m 

On the lattice, the Wilson term is not invariant under the rotation in eq. (J7J) and 
the twisted mass operator takes the form 

Am = £>w + [mo + i/^fCr 5 T 3 ] (10) 

where denotes the Wilson operator without the mass term, i.e. 

D°r = D w [m = 0] . (11) 

It can be shown that the twisted mass action leads to an 0(a)-improvement 
when the angle uj = tt/2. It goes beyond the scope of this paper to provide the 
arguments for this remarkable result the derivation of which can be found in ref. jH]. 
We only would like to remark that in order to obtain this value of the twist angle is 
equivalent to tune the bare fermion mass parameter m to a critical value m™ 1 *. 

One very important aspect of twisted mass fermions is that there is a particular 
definition of m" 1 ' from the vanishing of the PCAC fermion mass that does not only 
lead to an 0(a)-improvement j2H 123 UHl El EI but also substantially reduces cut-off 
effects that appear in 0(a 2 ) as has been demonstrated in [2S1IIH]- In this paper, 
we will use a value of m™ 1 that was obtained in the pure Wilson fermion theory 
without twisted mass term from the PCAC relation. The value of m™' was then 
tuned in such a way that the corresponding PCAC fermion mass vanishes. After 
having determined mjj nt in this way, we varied the twisted mass parameter fif to 
realize the fermion and pion masses we are interested in. 
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2.4 Overlap fermions 

A lattice Dirac operator -Dgw that satisfies the Ginsparg- Wilson relation |2Zj 

£>gw75 + 75^gw = 2aD GW -y 5 RD GW , (12) 

where R is a local term, leads to an action that has an exact (lattice) chiral symmetry 
eliminating thus automatically 0(a) cut-off effects. The realization of an operator 
-Dgw that we use here is the overlap fermion, which is characterized by the overlap 
Dirac operator [TU]. For R n , m = %^ it takes the form 

D ov = (l-^f-)D^+m ov , (13) 

D$} = ^{l + Do/y^} (14) 

where Dq is the so-called overlap kernel operator and a = a/ p. For the kernel Do, 
there is a large choice. In the following, however, we will only use the hypercube 
operator, i.e. D = D hyp with fixed A = 0.5, i.e. setting mh yp = — 1- This 
corresponds to setting the parameter p = 1 in eq. (|14j) and guarantees the locality 
of the overlap operator [2*%] . 

We remark that we have realized the square root operator 1/ J DJD q in eq. (J 14)) 
by Chebyshev polynomials P n e with degree n and a lower bound e. This Chebyshev 
polynomial shows an exponential convergence rate in the interval [e, 1]. Setting e 
to the lowest eigenvalue of the overlap kernel operator DqDq which is normalized 
to one, we always have chosen the degree n such that we reach machine precision 

for the evaluation of 1/v/ D\Dq. We used eigenvalue deflation and projected out a 
number of low-lying eigenvalues. In this case, e was set to the lowest non-projected 
eigenvalue of D\Dq. 

3 Observables and simulations 

In this section we give the operators that we have used to determine the physical 
observables we have computed and describe the numerical simulations. In particular, 
we discuss some of our attempts to perform dynamical overlap simulations avoiding 
problems with very low-lying eigenvalues of the overlap kernel operator. 

3.1 Observables 
2-point functions 

Generally, the bi-linear operators are given by 

T n T = ^ n rrVn, (15) 
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where T stands for certain combinations of the Pauli matrices a acting in Dirac- 
space, while the Pauli matrices r act in flavour space. We will consider the following 
operators, 



pa 
<Sn 



Aim 

yfia 
v n 



ifj n a 3 



(Dirac) (flavour) i 
'a M 



(Dirac) ^(Dirac) ^(flavour) 



r. 



(Dirac) ^(flavour) 



(16) 
(17) 
(18) 
(19) 



with V the pseudo scalar, S the scalar, A the axial-vector and V the vector operator. 

For the twisted mass fermions, the operators take a modified form as can be 
obtained from the field transformation according to eq. (J7J), leading to 



T>a(ph) 
' n 

A^a(ph) 
yp,a(ph) 



V, 



a(tb) 



3(tb) 



COS UJ V'i 
COS UJ An 



sin u; 6A 
f sin u e a bV 



ftb(tb) 



/in 



fi3(tb) 



•.'ifiaitb) 

cos uJVn + smuj e a bA 



fib(tb) 



u3(ph) 



if a 


= 1,2 


if a 


= 3 


if a 


= 1,2 


if a 


= 3 


if a 


= 1,2 


if a 


= 3 



(20) 
(21) 
(22) 



Here "ph" denotes the physical and "tb" the twisted basis. Note that in the special 
case uj = 7r/2 for a = 1,2, V and A just interchange their role while V remains 
invariant. In particular, the scalar operator is given by the 3rd component of the 
pseudo scalar operator. 

For the Wilson, hypercube and Wilson twisted mass fermions we computed the 2- 
point correlation functions by standard techniques using a conjugate gradient solver. 
For overlap fermions we calculated the correlators both from using a conjugate 
gradient solver and from eigenvectors <pi and eigenvalues A» of the overlap operator 
in eq. (|T3*j) . A generic non-singlet correlator C(n,m) is then expressed in terms of 
the fermion propagator S(n,m) and a suitable Dirac structure T, corresponding to 
the operators listed above, as 



C(n,m) = tr[rS(n,m)rS(m,n)} 



(0] Q (^)r a ,0f(n))(0r(m)r 75 05(m)) 



(23) 



PCAC fermion mass 

The bare mass parameter for Wilson (m ) and hypercube {m^y P ) fermions re- 
ceives an additive mass renormalization, necessitating the determination of a critical 
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fermion mass m cr i t which we computed by the vanishing of the PCAC fermion mass, 
extracted from the PCAC relation 



{d^A"(n)d) = 2m(V(n)d). (24) 

Choosing O = V and projecting to zero momentum (denoting by A and V the 
currents summed over space), we arrive at 

m f (t) = (25) 

For the currents that appear in eq. (J23J) one may use the local currents of eqs. (fTHl - 
[TT)|h An alternative is to employ also conserved currents. A very general way to 
derive the conserved currents is given in [2H] . The vector and axial currents are then 
given by 

V£ a = jf + J» L (26) 
A? = Jf-J£ L (27) 



where 



^ = ^(1-^)^(1 + ^ (28) 
J^ = ^(l + a 3 )/C^(l-^ (29) 



and 



a=0 



i ^ n,fi 



lM = e^Un^ (30) 



where Dill^) denotes the lattice Dirac operator used. This method of constructing 
the conserved currents is very useful when complicated lattice Dirac operators are 
considered since their construction from the current conservation condition can be 
very cumbersome. We followed this prescription to compute the conserved currents, 
only for the overlap fermions we used the local point currents. We remark that for 
the twisted mass and the overlap case in principle also the bare fermion masses, 
fif and m ov can be used. However, since this could lead to very different lattice 
artefacts, we employed also in these cases the fermion mass derived from the PCAC 
relation for all physical results presented in the following in order to be able to 
directly compare the different fermion actions. 

Scalar condensate 

Besides the 2-point functions which will provide the pseudo scalar mass, we 
considered also the scalar condensate, E = (ipip). A first method to compute E is 
by calculating TiD^ 1 using Gaussian noise sources. We will denote the so computed 
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values of £ as ^direct- This quantity develops in the case of Wilson and hypercube 
fermions a divergent piece oc 1/a. A second way, which avoids the appearance of 
the divergent piece from the beginning, is to use the integrated axial Ward identity 
leading to a "subtracted" scalar condensate, £ su b jSDj- We remark that in the case 
of overlap fermions we computed £ from the improved scalar operator (^(1 — ^D)ijj) 
which we evaluated from the the eigenvalues of the overlap operator. More precise 
definitions and a further discussion will be provided in section 4. 



3.2 Simulations 



In our work, we have used for the Wilson, the hypercube and the twisted mass 
fermion action a standard Hybrid Monte Carlo algorithm (HMC) [HI]. In the case 
of overlap fermions a straightforward implementation of the HMC algorithm is not 
suitable since the variation of the Hamiltonian with respect to the gauge field can 
lead to very large forces proportional to (-DJ-Dq)" 3 / 2 when the overlap kernel operator 
Dq develops exceptionally small eigenvalues. See refs. [H21 EHl HH] for variants of the 
HMC algorithm that circumvent this difficulty 

For the overlap fermion results we used configurations that were generated in the 
pure gauge theory only and performed a reweighting with the overlap fermion de- 
terminant to obtain physical observables for dynamical Nf = 2 flavours of fermions. 
However, we also tried several possibilities to avoid the problem with low-lying eigen- 
values of the overlap kernel operator D by replacing D ov by some operator D^ pprox 
that is a good approximation to D ov but which is safe against these low-lying eigen- 
values. The simulation can then be made exact again by adding a correction step 
employing the ratio D ov /D^ ppmx . 

The general idea is to write (we use for simplicity only a single operator here, in 
practice one would have to use the operator D^D, of course) 



det A- 



detD^ pmx ■ det 



D r 



r-\approx 

Una 



det D™ • ^dct • 



(31) 



While detD'o Ppmx would be used in the HMC algorithm, the remaining determinant 
ratio Rdct could be implemented either in an additional accept/reject step or it could 
be included as a reweighting factor in the computation of a given observable. 

A crucial question in such an approach is, whether an operator D^ pprox can 
be found such that the fluctuations in i?det are small enough to obtain statistically 
significant results. We decided therefore to test this idea by computing R dct stochas- 
tically using n — 10 Gaussian noise vectors on a number of gauge field configurations 
generated in the pure gauge theory at (3 = 3 on a 16 2 lattice. In the following we 
will describe the results of these tests for three choices of D^ pprox . 



Case of D* Ppiox 



D 



hyp 



As a first trial, we used the hypercube operator -Dhyp as an approximation to D ov . 
This choice has been motivated by the fact that the operator -Dhyp is constructed to 
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be approximating the Ginsparg- Wilson relation, resulting in a very similiar eigen- 
value spectrum for both operators ^3] . Note also that we have used the hypercube 
operator itself as an overlap kernel operator Dq. Let us remark that we re- write the 
massive overlap operator as D$ + m ov /(l — m ov /2) in order to better match the 
eigenvalue spectra of both operators. 

We computed R^t as a function of the hypercube bare fermion mass m^ yp in a 
range —0.25 < mh yp < 0.25. However, we found that the fluctuations of i?det were 
extremely large for all values of mh yp we have tested. Thus we had to conclude 
that -Dhyp cannot serve as an infrared safe, approximate operator for the overlap 
simulations when R det is used in a stochastic correction step. We cannot exclude, 
of course, that by using an improved hypercube operator, e.g. by adding the clover 
term or performing smearing of the link variables, the situation could be improved. 
However, given the negative findings of our investigation, we did not pursue this 
direction. What we tried instead, is to use -Dhyp as the guidance Hamiltonian in the 
molecular dynamics part of the Hybrid Monte Carlo algorithm while keeping D ov as 
the exact overlap operator for the accept/reject Hamiltonian. Amazingly, despite 
the negative results described here for i?det, we found that this led to reasonable ac- 
ceptance rates, at least for not too large systems, see ref. jUS] for a further discussion 
of this point. 

Modified Chebyshev polynomial 

As a second attempt, we tried to modify the range of the Chebyshev polynomial 
employed in the construction of the overlap operator. There, the square root of the 
kernel is computed by a Chebyshev polynomial P n)6 of degree n in the interval [e, 1], 



We have chosen the approximation accuracy in eq. (|32jl to be very high, compatible 
with machine (64-bit) precision. The value of e is chosen to correspond to the 
lowest eigenvalue of D^Dq. In principle, it is possible to use instead of the "exact" 
polynomial with parameters n and e a modified polynomial with parameters n < n 
and e > e. Clearly, if e could be chosen to be well above the smallest eigenvalues 
of DlD then using an overlap operator with a polynomial g in the molecular 
dynamics part would lead to an infrared safe simulation. 

As a start situation, we had chosen e = e and reduced only the degree of the 
polynomial successively. We found that it is indeed possible to reduce n substantially 
to n ~ 1/3-n without having large fluctuations in the determinant ratio -Rdet- This is 
a very positive outcome since it gives rise to a substantial acceleration of simulations 
with the overlap operator. However, when we tried to change the value of e only 
slightly, we observed immediately large fluctuations in -Rdet- For e = 1.1 • e the 
fluctuations in R^ ct became already so large that there is no chance to obtain a 




(32) 
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statistically significant result in such a setup. 

Explicit infrared regularization 

As a third (and last) attempt, we tried to use an approximate overlap operator 
that has an infrared regulator built in. In particular, we studied the situation where 
we use a modified sign function 



The (optimistic) expectation here is that by choosing 5 large enough, without having 
too large fluctuations in the determinant ratio Rdet, the overlap operator with such 
a modified sign function would be infrared safe for the HMC simulation. 

Such an optimistic expectation is not completely unfounded when one inspects 
the eigenvalues of the exact overlap operator D ov and the sign function modi- 
fied one, D^ pro *. Indeed, we show in fig. ^ the eigenvalues when choosing 5 = 
0, 0.0001, 0.001, 0.01 for a typical configuration. 

As can be seen, all the eigenvalues lie on the expected circle and a difference for 
various choices of 6 is not visible. In fig. |2]we show the difference of the eigenvalues 
between the spectra using 5 = and 5^0. In building the difference of the 
eigenvalues we used an ordering of the eigenvalues with respect to their real part. 
The scale in the plot is chosen for the situation of 5 = 0.01. We have rescaled the 
difference in the eigenvalues by a factor of 30 for 5 = 0.0001 and 7 for 5 = 0.001 in 
order to plot all difference spectra in one common plot. It appears that the difference 
spectra do not build a perfect circle shape but are rather lemon shaped. Nevertheless 
the lemon distortion of the circle happens at a rather small scale of O(10~ 2 ) for the 
case of 5 = 0.01 while for smaller values of S the difference becomes even smaller. 
Note that the smallest eigenvalue is A min (.DoA)) = 0.24 for this configuration. 

The smallness of the distortions in the spectrum as observed in fig. |21 appears to 
be quite promising for this choice of an approximate overlap operator. However, the 
effects on the fluctuations are considerable as can be seen in fig. EJ For a value of 
5 = 0.0001 (circles) the determinant ratio i?det ~ 1 and has very small fluctuations 
until the degree of the polynomial is lowered to n = 8 which is a factor of three 
smaller than the original degree of the polynomial. However, already for 5 = 0.001 
(upward triangles) i?det starts to develop some fluctuations and finally, for 6 = 0.01 
(downward triangles) the fluctuations are becoming so large that realistic simulations 
with such a value of 5 cannot be performed. 

As a conclusion, we found that none of the above described attempts led to 
a satisfactory solution for an infrared safe dynamical overlap simulation. As said 
above, we therefore have performed the overlap simulations finally by generating 
pure gauge configuration and using a reweighting with the determinant, computed 




(33) 
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1.5 




_]__5 1 1 1 1 1 1 1 

-0.5 0.5 1 1.5 2 2.5 

Real Z 

Figure 1: The eigenvalues of the exact overlap operator and a modified overlap oper- 
ator when the sign function is changed according to eq. (jHSJ)- The plot contains the 
eigenvalues for 5 = (the original overlap operator) and for 5 = 0.0001, 0.001, 0.01. 
For all these cases, the eigenvalues are so close to each other that they cannot be 
distinguished in the graph. 
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0.015 




-0.015 ' 1 1 1 1 1 1 

-0.015 -0.01 -0.005 0.005 0.01 0.015 

Real Z 

Figure 2: The difference spectra of the overlap operator with 5 = and the operators 
with 5 = 0.01 (most outer spectrum), 5 = 0.001 (middle spectrum) and 5 = 0.0001 
(most inner spectrum). The difference spectra for 5 = 0.001 and 5 = 0.0001 have 
been scaled with a factor 7 and 30, respectively, in order to be able to include them 
in the plot. 
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10 




Figure 3: The determinant ratio R^ et as a function of the degree n of the Chebyshev 
polynomial. We plot the cases for 5 = 0.0001 (circles), 5 = 0.001 (upward triangles) 
and S = 0.01 (downward triangles). 



exactly from the eigenvalues. 

Summary of simulation setup 

In the following table we shortly summarize which technique we have used for 
the different lattice fermions for the simulations and for the observables. 





Wilson / tm / hypercube 


overlap 


simulations 


standard HMC 


determinant re weighting 


operators 


conjugate gradient solver 


eigenvalues 

& eigenvectors/ 
conjugate gradient solver 


currents 


local / conserved 


local 


pseudo scalar correlator 


(VV) 


(VV - SS) 


scalar condensate 


Ward-Takahashi-identity 


(V(l - f D)V) 



Let us end this section by a small technical remark. For the computations of 
the eigenvalues and eigenvectors we used the LAPACK [3B] routine. In the course 
of our work we found that this routine does not compute correctly the eigenvectors 
in the case of degenerate zero modes. It rather gives linear combinations of the 
exact solutions which led to a problem in the computation of the pseudo scalar 
correlator using the operator V if configurations with topological charge \Q\ > 1 
are considered. The effect shows up in such a way that the pseudo scalar correlator 
eq. (fTH|) develops a plateau like behaviour for values of Euclidean time t close to the 
middle of the lattice. This behaviour made it practically impossible to extract the 
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pseudo scalar mass and the PCAC fermion mass. By re-diagonalizing the zero mode 
sector we could verify that this effect goes away. 

As a simple solution for our simulations we decided to always take the correlator 
(V(t)V(0) — S(t)S(0)) in which this zero mode contribution is cancelled out |57| . 
The only exception is the case z = 0.8, where the pseudo scalar masses are so 
high that it became very difficult to disentangle them from the scalar masses when 
the (iS(t)iS(O)) was subtracted of. Therefore we used for this large value of z the 
conjugate gradient solver to compute the correlators. 



4 Results 



It is the main aim of this paper to check the scaling behaviour of the fermion actions 
described in section 2. To this end, we have chosen the simulation parameters such 
that we fix the scaling variable 



i m f/g) 



2/3 



m 



2/3 



(34) 



We performed simulations for a wide range of j3- values, 0.1 < < 6. For Wilson, 
hypercube and twisted mass fermions our lattices were mainly of size V = L 2 = 32 2 . 
Only at /3 = 5, 6 and z = 0.2 we went up to V = 48 2 lattices. For the overlap 
fermion simulations we mainly used V = 20 2 and, for z = 0.2, V = 24 2 lattices. 
Only at (3 — 5 and z = 0.2 we used a V = 28 2 lattice. It is our experience that for 
larger lattices the determinant reweighting technique used for the overlap fermion 
simulations are not practical since they lead to large fluctuations which spoil the 
signal to noise ratio. We finally mention that all our error analyzes are based on 
the method described in ref. [HE] including thus the autocorrelation times in the 
computation of the errors. 

As a prime quantity we will test the scaling behaviour of m n \/~/3, performing 
finally a continuum limit for all actions used. In the continuum, there are two 
approximate calculations for m™ nt /e for the massive Schwinger model. 

The first is performed by Smilga [39J who finds for strong coupling and small 
fermion mass 



cont 



g5/6 r7 /3/ r(3/4) y/3 r(i/6) /my 
vr(i/4)/ r(2/3)V e y 



x 2/3 



fm\ 2 / 3 
= 2-008 ( 7 ) • 

For large masses, Gattringer in ref. [10] finds from a semi-classical analysis 



m 



cont 
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(35) 
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Wilson 
"^crit 


hyp 
- m crit 


1.0 


0.3204(7) 


0.335(1) 


2.0 


0.1968(9) 


0.203(1) 


3.0 


0.1351(2) 


0.1392(4) 


4.0 


0.1033(1) 


0.1050(2) 


5.0 


0.0840(1) 


0.0856(1) 


6.0 


0.0719(1) 


0.0727(1) 



Table 1: Critical hopping parameters at the values of (3 used in our simulations for 
the Wilson, m™ Bon , and hypercube, m^, fermions. 



Note that we give these results in terms of the dimensionful quantities, coupling e, 
fermion mass m and pseudo scalar mass m™ nt . Both expression are approximate 
computations and it is interesting to compare these against our non-perturbative 
calculations. 

For completeness, we also give here analytical expressions for the scalar conden- 
sate as available in the literature. The first is again from ref. |39j . 

S cont /m\V3 

« 0.388 ( — ) . 39 

e V e / 

A second expression is derived in ref. |4*T| . 

— = T- f A • ( 4 °) 
e 4tt (m/e) 



Critical fermion mass 

For Wilson fermions and for hypercube fermions, the fermion mass receives an ad- 
ditive renormalization. Therefore, the critical values of the bare fermion mass, where 
e.g. the PCAC fermion mass vanishes has to be determined by non-perturbative sim- 
ulations. We computed mpcAC; extracted from eq. ()25|) using the conserved currents 
of eqs. ()26I27|) as a function of the bare fermion masses, mo for Wilson and mh yp for 
hypercube fermions. We determined at each value of /3 those values of mo and mh yp 
where mpcAC — 0. The value of these bare fermion masses can be found in tabled 
Note that the actual critical mass values used for the twisted mass simulations in 
order to realize full twist, though compatible, differ slightly from the ones in tabled 
since they were obtained with less statistics. 

Fermion mass dependence of the pion mass 

As described above, in the scaling analysis we will fix the scaling variable z = 
(m/v^) 2 ^ 3 - To this end, we first explored the dependence of the pion mass as a 
function of the PCAC fermion mass for fixed values of (3. We give one example for 
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0.2 0.4 0.6 0.8 1 

(m/^) 2/3 

Figure 4: Pseudo scalar mass m^y/P versus z = (mf\f]3) 2 ^ 3 for three values of (5. 
The solid lines are the analytical predictions of eqs. ()35I37|) . Crosses represent our 
data for [3=1, triangles for /3 = 0.5 and diamonds for (3 = 0.1. 

the case of hypercube fermions in fig. fusing the conserved currents of eqs. (|26I27|) to 
compute the PCAC fermion mass. We also plot in the graph the analytic predictions 
of eqs. (|35I37|) for the fermion mass dependence of the pion mass. 

The data in fig. H are for (3 = 0.1, (3 = 0.5 and (3 = 1.0. Clearly, for j3 = 0.1 
the data do not follow the theoretical expectation while for the other values of /3 
there seems to be some agreement with the analytical formulae. However, in order 
to make more definite statements, a closer look to the scaling behaviour is clearly 
needed. Anyhow, the main purpose of fig. 0] is to show that we have collected data 
that are concentrated around the values of z = 0.2, z = 0.4 and z = 0.8 that we are 
interested in for our final scaling analysis. When the values of z did not coincide 
directly with the desired ones, we performed an only very small linear interpolation 
of our data to achieve the exact values of z. 

Finite size effects 

One source of a systematic error in the determination of are possible finite 
size effects. For two dimensions the asymptotic finite size corrections for the pseudo 
scalar mass were computed in ref. [12] and studied numerically in the case of the 
Schwinger model in ref. [3]. A very good agreement to the theoretical prediction 
was found and it was observed that the variable m^L needs to be surprisingly large 
to suppress finite size corrections. 

Taking this as a warning, we followed the procedure of ref. [3] and performed a 
finite size correction of our values for the pseudo scalar masses when necessary. We 
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T 




Figure 5: We show m- K \[ft as a function of L at the example of Wilson fermions 
and three values of (3 (downward triangles: (3 = 4, circles: f3 — 5, upward triangles: 
(3 = 6). The fit curves, represented by solid lines, are obtained by performing 
a common fit to eq. (j41j) to our data for z = 0.2 and z = 0.4 and for Wilson, 
hypercube and twisted mass fermions, choosing always (3 > 4. 

used the analytical formula 



m 7r (L) = rrC + A^^e- L ^ (41) 

where denotes the pseudo scalar mass in the infinite volume limit. Rescaling 
eq. (j4~Tj) by y//3 this formula can be applied to all our finite volume data at various, 
large enough values of (3 to keep the effects of the lattice spacing small. Inspecting 
our data, we decided to perform a global fit to our simulation data obtained with 
Wilson, hypercube and Wilson twisted mass fermions at z — 0.2 and z = 0.4 for 
(3 > 4. As a result we found a "universal" constant A = 12.4(6) which is compatible 
with the value computed in ref. [3]. In fig. El we give an example for the resulting fit 
for a subset of our data in the case of Wilson fermions. 

We then examined all pseudo scalar masses that we will use in the detailed 
scaling analysis described below and used eq. (|41|) to analytically correct for finite 
size effects. We give in tables |2J El and 0] the infinite volume values for the pseudo 
scalar mass as obtained from this finite size correction. Note that in most of the 
cases, this correction is negligibly small and stays anyhow on the percent level. 

4.1 Results for the pseudo scalar mass 

In tables El El and El we give our results for the pseudo scalar masses at fixed scaling 
variable z of eq. (J34|) . where z is fixed either by using local or conserved currents to 
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Wilson fermion results 



z 








,™cons 



L 


local m 

7T,00 V f 


cons m 

7T,00 V f 


£ c T\//3 

sub V f 


0.2 


1.0 


-0.231367 


-0.230193 


32 


0.388(5) 


0.391(3) 


0.277(7) 


0.2 


2.0 


-0.132316 


-0.131842 


32 


0.402(3) 


0.406(2) 


0.252(6) 


0.2 


3.0 


-0 082626 


-0 081935 

V 7 • V / V 7 -L '.7 ' J .J 


32 


0.404(4) 


408(2) 


232(7) 


0.2 


4.0 


-0.057374 


-0.056979 


32 


0.406(2) 


0.408(2) 


0.231(8) 


0.2 


5.0 


-0.043199 


-0.043010 


48 


0.403(2) 


0.404(1) 


0.219(6) 


0.2 


6.0 


-0.034249 


-0.034009 


48 


0.404(2) 


0.406(2) 


0.225(4) 


0.4 


1.0 


-0.075308 


-0.059286 


32 


0.75(1) 


0.779(5) 


0.339(6) 


0.4 


2.0 


-0.017495 


-0.008352 


32 


0.783(8) 


0.813(6) 


0.339(3) 


0.4 


3.0 


0.014505 


0.019789 


32 


0.799(7) 


0.819(3) 


0.334(4) 


0.4 


4.0 


0.026358 


0.031078 


32 


0.80(1) 


0.820(3) 


0.340(3) 


0.4 


5.0 


0.032526 


0.035391 


32 


0.808(4) 


0.824(2) 


0.345(2) 


0.4 


6.0 


0.036146 


0.037693 


32 


0.811(3) 


0.820(2) 


0.345(2) 


0.8 


1.0 


0.352443 


0.537801 


32 


1.36(1) 


1.560(5) 


0.329(1) 


0.8 


2.0 


0.314017 


0.403487 


32 


1.52(1) 


1.690(5) 


0.400(1) 


0.8 


3.0 


0.298677 


0.350958 


32 


1.60(1) 


1.729(4) 


0.438(1) 


0.8 


4.0 


0.276568 


0.314133 


32 


1.63(1) 


1.745(5) 


0.469(1) 


0.8 


5.0 


0.255009 


0.285628 


32 


1.65(1) 


1.758(3) 


0.494(1) 


0.8 


6.0 


0.239625 


0.262974 


32 


1.686(6) 


1.770(3) 


0.518(1) 



Table 2: Results for the Wilson fermion simulations. The scaling variable (mji/^) 
is fixed by using either the local or the conserved currents to determine the PCAC 
fermion mass resulting in the two different values of the pseudo scalar mass. Typical 
statistics of the runs were between 4000 and 10000 configurations. 
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Hypercube fermion results 



A/ 


ft 




,™coiis 


r 




m coas ,fB 


Vcons m 
^sub VP 


n 9 


1 n 




-0 931 D44 


39 


3QH 1 






n 9 


9 n 


-0 1 37038 

U .10 1 KJOO 


-D 1 333RQ 


39 


D 4(Y)('3l 


D 41fif3l 

U.tIU I O 1 




n 9 

U. Zi 


3 n 




-fl 0831 QQ 


39 


U.t:UU 1 I J 


D 41 5f9l 




0.2 


4.0 


-0.059364 


-0.058095 


32 


0.399(9) 


0.416(2) 


0.240(5) 


0.2 


5.0 


-0.044387 


-0.043885 


48 


0.408(2) 


0.410(2) 


0.246(8) 


0.2 


6.0 


-0.035528 


-0.034893 


32 


0.413(7) 


0.411(7) 


0.247(4) 


0.4 


1.0 


-0.083223 


-0.045098 


32 


0.78(2) 


0.844(9) 


0.375(2) 


0.4 


2.0 


-0.025253 


-0.004337 


32 


0.78(1) 


0.846(4) 


0.381(2) 


0.4 


3.0 


0.012877 


0.021147 


32 


0.806(6) 


0.836(3) 


0.382(2) 


0.4 


4.0 


0.025787 


0.031202 


32 


0.809(4) 


0.834(4) 


0.385(2) 


0.4 


5.0 


0.031525 


0.035539 


32 


0.812(3) 


0.831(2) 


0.392(3) 


0.4 


6.0 


0.034082 


0.037392 


32 


0.816(3) 


0.835(2) 


0.392(3) 


0.8 


1.0 


0.311570 


0.613589 


32 


1.39(2) 


1.742(6) 


0.366(1) 


0.8 


2.0 


0.309020 


0.433230 


32 


1.55(2) 


1.792(6) 


0.466(1) 


0.8 


3.0 


0.293918 


0.365237 


32 


1.61(1) 


1.793(4) 


0.517(1) 


0.8 


4.0 


0.272609 


0.321868 


32 


1.65(1) 


1.792(5) 


0.557(1) 


0.8 


5.0 


0.255540 


0.290832 


32 


1.672(7) 


1.792(4) 


0.588(1) 


0.8 


6.0 


0.237752 


0.266771 


32 


1.68(1) 


1.796(5) 


0.616(1) 



Table 3: Results for the hypercube fermion simulations. The scaling variable 
(m/v7?) 2 ^ 3 is fixed by using either the local or the conserved currents to determine 
the PCAC fermion mass resulting in the two different values of the pseudo scalar 
mass. Typical statistics of the runs were between 4000 and 10000 configurations. 
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Twisted mass fermion results 


z 


P 




™,cons 


L 




"CSV? 


^sub V P 


0.2 


1.0 


- 


0.071591 


32 


- 


0.434(2) 


0.303(7) 


0.2 


2.0 


0.050499 


0.054990 


32 


0.408(8) 


0.428(2) 


0.271(10) 


0.2 


3.0 


0.045007 


0.046886 


32 


0.411(2) 


0.422(2) 


0.249(6) 


0.2 


4.0 




0.041533 


32 




0.421(1) 


0.226(6) 


0.2 


5.0 


0.036965 


0.037643 


48 


0.409(2) 


0.414(2) 


0.237(5) 


0.2 


6.0 




0.034651 


48 




0.414(2) 


0.233(7) 


0.4 


1.0 


0.177303 


0.208213 


32 


0.78(1) 


0.874(3) 


0.446(7) 


0.4 


2.0 


0.144653 


0.157190 


32 


0.81(1) 


0.867(3) 


0.412(4) 


0.4 


3.0 


0.125186 


0.133453 


32 


0.819(8) 


0.851(2) 


0.394(4) 


0.4 


4.0 


0.113557 


0.117956 


32 


0.824(4) 


0.846(2) 


0.385(5) 


0.4 


5.0 


0.102819 


0.106827 


32 


0.820(4) 


0.844(2) 


0.388(4) 


0.4 


6.0 


0.095373 


0.098326 


32 


0.821(3) 


0.842(2) 


0.393(5) 


0.8 


1.0 


0.471381 


0.654028 


32 


1.47(2) 


1.799(5) 


0.571(2) 


0.8 


2.0 


0.376485 


0.459309 


32 


1.60(1) 


1.856(3) 


0.628(2) 


0.8 


3.0 


0.334417 


0.384597 


32 


1.66(2) 


1.838(3) 


0.638(2) 


0.8 


4.0 


0.306509 


0.337947 


32 


1.695(6) 


1.835(3) 


0.659(2) 


0.8 


5.0 


0.282663 


0.305270 


32 


1.719(6) 


1.829(3) 


0.671(2) 


0.8 


6.0 


0.262653 


0.280421 


32 


1.731(6) 


1.833(3) 


0.683(2) 



Table 4: Results for the Wilson maximally twisted mass fermion simulations. The 
scaling variable {mf^/]3) 2 ^ is fixed by using either the local or the conserved currents 
to determine the PCAC fermion mass resulting in the two different values of the 
pseudo scalar mass. Typical statistics of the runs were between 4000 and 10000 
configurations. 
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Overlap fermion results 



p 


m ov 


(movvW 73 


(m^'v^) 273 


L 






3.0 


0.0464758 


0.1864 


0.194(15) 


24 


0.37(4) 


0.176(4) 


4.0 


0.0447214 


0.2000 


0.200(8) 


24 


0.42(4) 


0.177(4) 


5.0 


0.0420000 


0.2066 


0.20(3) 


28 


0.46(4) 


0.138(2) 


1.0 


0.2529822 


0.4000 


0.4061(14) 


20 


0.807(16) 


0.1950(16) 


2.0 


0.1788854 


0.4000 


0.4069(19) 


20 


0.824(14) 


0.2340(12) 


3.0 


0.1460595 


0.4000 


0.399(2) 


20 


0.851(12) 


0.2518(7) 


4.0 


0.1264911 


0.4000 


0.383(2) 


20 


0.833(10) 


0.2618(7) 


5.0 


0.1142685 


0.4027 


0.396(5) 


24 


0.837(10) 


0.2415(8) 


1.0 


0.5724334 


0.6894 


0.797(2) 


20 


1.433(10) 


0.2363(3) 


2.0 


0.4232769 


0.7103 


0.8047(5) 


20 


1.568(4) 


0.2905(3) 


3.0 


0.3553500 


0.7236 


0.7934(2) 


20 


1.604(3) 


0.3270(3) 


4.0 


0.3152474 


0.7353 


0.7916(1) 


20 


1.638(2) 


0.3549(3) 


5.0 


0.2916839 


0.7521 


0.8010(2) 


20 


1.682(2) 


0.3811(6) 



Table 5: Results for the overlap fermion simulations. The scaling variable 
(m/v7?) 2 ^ 3 is fixed by using the local currents to determine the PCAC fermion 
mass. We give for comparison also the values of the scaling variable z when the 
bare overlap fermion mass m ov is used. Typical statistics of the runs were 4000 
configurations. 
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Figure 6: We show m^^ffi as a function of (3 at fixed value of z = 0.4 as determined 
from the PCAC fermion mass employing conserved currents. We represent results 
from Wilson fermions by circles, from twisted mass fermions by downward and 
for hypercube fermions by upward triangles. The arrows represent the theoretical 
prediction, the upper arrow is from eq. (|37|) the lower from eq. (|35|). 



extract the PCAC fermion mass. While in the tables we give only the results for 
(3 > 1, we have performed simulations also for smaller values of (3 for all actions, 
except overlap fermions. We show the result for m n y/j3 as a function of (3 in fig. El 
The graph is done for the example of z = 0.4 where we used the conserved currents to 
determine the PCAC fermion mass. Clearly, large lattice cut-off effects are observed 
when small values of (3 < 1 are taken. From the figure it is obvious that for an 
asymptotic scaling analysis values of (3 > 1 should be taken. 

In figs. |H1 and El we show the results of our scaling test for fixed z = 0.8, 
z = 0.4 and z = 0.2 respectively. For the values of (3 — 1 and (3 = 2 at z = 0.2, the 
determinant used for reweighting the overlap results induced large fluctuations such 
that no reliable values of physical observables could be extracted. For all figures 
the value of z was determined using local currents to compute the PCAC fermion 
mass. We show m n y/j3 as a function of 1/(3 oc a 2 . We first performed linear fits 
in 1/(3 independently for each lattice fermion used. Since these fits gave consistent 
continuum values for m v ^f(3 for all actions, see the example for z = 0.4 in ref. we 
finally performed constraint fits, demanding that all actions give the same continuum 
value for m n y/]3 at fixed value of z. We show these constraint fits in figs. [Zl [HI and El 
as the solid lines. The data for all kind of fermions, Wilson, hypercube, maximally 
twisted mass and overlap, are nicely consistent with a linear behaviour in a 2 . While 
this is expected for maximally twisted mass and overlap fermions, this outcome is 
somewhat surprising for hypercube and Wilson fermions. Note that our finding for 



22 




Figure 7: Scaling test of m n y/]3 as function of 1/(3 oc a 2 at fixed value of z = 0.8. 
The solid lines represent linear constraint fits in 1/(3 demanding the same continuum 
limit value for all kind of lattice fermions used. For Wilson, hypercube and twisted 
mass fermions we used data obtained at (3 > 3, while for overlap fermions data with 
(3 > 2 were used. 



Wilson fermions is, however, consistent with the results in ref. |43j . 

In general it is very difficult to decide which kind of lattice fermion shows the 
best scaling behaviour and we cannot draw a definite conclusion here. First of all, 
all kind of fermions show only small lattice artefacts. Second, if one would use the 
PCAC fermion mass from the conserved currents, the picture changes. We give an 
example for z = 0.4 in fig. El which reveals again an 0(a 2 ) scaling for the here 
used Wilson, maximally twisted mass and hypercube fermions. When compared to 
fig. |H1 where maximally twisted mass fermions show the smallest scaling violations, 
in fig. hypercube fermions seem to do better. 

Continuum comparison to theory 

As was shown above, all kind of fermions show a nice scaling behaviour that is lin- 
ear in a 2 and give a universal continuum limit. In fig. ^Jwe show m^^f^ as a function 
of z in the continuum and compare to the theoretical expectations of eq. (J33j) (lower 
line) and eq. (}3Tj) (upper line). As an inlay we plot the ratio R m ^ = m^ ata / m^ heor 
of our non-perturbatively obtained data extrapolated to the continuum and the two 
theoretical predictions. Triangles represent the data divided by the corresponding 
value computed from eq. (jBTj) while circles represent the data divided by the cor- 
responding value computed from eq. (|33j) . To the precision we could obtain in this 
work, the theoretical predictions do not describe the non-perturbatively obtained 
simulation data satisfactory at all values of z. Only for small values of the fermion 
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Figure 9: Scaling test of m n y/J] as function of 1/(3 oc a 2 at fixed value of z = 0.2. 
The solid lines are fits explained in fig. 
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Figure 10: Scaling of m^^fji as function of 1/(5 using conserved currents to compute 
the PCAC fermion mass to fix z = 0.4. We denote by circles Wilson, by downward 
triangles twisted mass and by upward triangles hypercube fermions. The solid lines 
represent linear fits in 1//3 for (3 > 3. The arrows represent the theoretical prediction, 
the upper arrow is from eq. (|37|) the lower from eq. (|35|) . 

mass (z = 0.2) there seems to be some agreement with eq. (|3*B]) while at z = 0.8 
eq. (|3T|) seems to hold only. 

4.2 Scalar Condensate 

Another physical quantity we will consider in this work is the scalar condensate 
E = (ipip), for which analytical predictions exist, see eq. ()39|) and eq. (|40[). A very 
simple way to calculate the scalar condensate is to compute 

S = ^^Tr,D- 1 (x,x) (42) 

X 

using a stochastic method. We denote the so computed values of E as E,-i irect . A 
severe drawback of this definition of E is that, at least in the case of Wilson and 
hypercube fermions, from the mixing with the identity operator a divergent piece 
ocl/a appears that needs to be subtracted non-periurbatively. 

In the case of the twisted mass fermions at full twist, it is possible to use the op- 
erator {2pa 3 T 3 ip} , i.e. the 3rd component of the pseudo scalar operator (see eq. (HEJ)), 
for the calculation of the scalar condensate. This operator does not mix with the 
identity operator and thus the divergent piece does not appear. 
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(m /x //?) 2/3 

Figure 11: Continuum behaviour of m n y/]3 as function of z. The lower curve rep- 
resents eq. (j35|) . the upper curve eq. (pTfjl . The inlay gives the ratio R m7V of our 
non-perturbatively obtained data and the two theoretical calculations. Triangles 
use eq. (137)) and circles use eq. (JH^J. 

Also for the overlap operator there is a definition S ov = (^(1 — ^D^ip) which 
subtracts of the divergent piece automatically. For overlap fermions, we calculated 
therefore directly the scalar condensate from the eigenvalues of D ov 

Sov^C^l-MO)/^), (43) 

i 1 

where \ denotes an eigenvalue of the massive overlap operator and Aj(0) denotes 
an eigenvalue at zero fermion mass, i.e. the eigenvalues of Dov . 

A second way, which avoids the appearance of the divergent piece from the 
beginning, is to compute a so-called subtracted scalar condensate S sub using the 
integrated axial Ward-Takahashi identity |3Uj . 

S sub = 2m PCAC ^(V X V Q ) . (44) 

In the twisted mass case, this corresponds to the integrated PCVC relation 

Ssub = (^e>"3T3^)sub = ^mpcvc^CPxVo) ■ (45) 

X 

Comparison of Sd ircct and E su b 

To just illustrate the effect of using a direct and a subtracted definition of the 
scalar condensate, we plot ^direct and S su b for Wilson and maximally twisted mass 
fermions in fig. 
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Figure 12: The condensate from the direct Street (°P en symbols) and the subtracted 
E su b (full symbols) determinations. Error bars are within the size of the symbols. 
Upper graph: S direct and S sub for Wilson fermion. We indicate the values of (3 where 
^direct ^ s obtained. Lower graph: S direct and S sub for twisted mass fermions. In this 
case, the results from both condensates are consistent up to scaling violations. 

In two dimensions, there is in general a relation between Street an< ^ ^sub given 

by 

S su b = S direct — c / a . (46) 

The coefficient cq multiplying the 1/a divergence is to be determined non-perturba- 
tively. In the upper panel of fig. EH the 1/a dependence is clearly visible for E direct , 
as for increasing values of (3, the values of Edirect increase accordingly. It is clear 
from the figure that an extraction of a physical value for the scalar condensate will 
be very difficult since the term Co/a dominates the signal. 

According to [20], Co comes from the explicit breaking term of chiral symmetry 
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for Wilson fermions. Therefore, one can expect that for twisted mass fermions at 
maximal twist this term is absent and Street behaves like X su t>. This is shown in the 
lower panel of fig. In the case of twisted mass fermions both Street an d ^sub are 
comparable (modulo scaling violations). In fact, we find a tendency that the ratio 
S su b/Sdircct approaches one when (3 is increased. We finally remark that also for the 
improved definition of the scalar condensate of eq. (J43|) in case of overlap fermions 
the divergence term is automatically subtracted of. As a result of this discussion 
we will calculate the scalar condensate from the subtracted scalar condensate in the 
case of Wilson, hypercube and maximally twisted mass fermions while for overlap 
fermions we will use the direct calculation with the improved definition. 

Scaling of the scalar condensate 

Let us now turn to the results of the scaling behaviour of the scalar condensate 
for the various fermion actions. We show the results in fig. ED S su b is calculated 
for Wilson, hypercube and twisted mass fermion, and Street is done for overlap 
fermions, at z = (m /v //3) 2/3 = 0.2, 0.4 and 0.8. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 




Figure 13: Scaling of the scalar condensate as a function of 1/(3 at z — 0.2 (circles), 
z = 0.4 (upper triangles) and z = 0.8 (downward triangle). The solid lines represent 
the fits using eq. (J4"%j) for (3 > 2. 

The figures show a strong dependence of the scalar condensate on in par- 
ticular for the heavy fermion region. This behaviour can be explained by the fact 
that in two dimensions, the scalar condensate develops a logarithmic divergence in 
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(3 when the lattice spacing is sent to zero *. This is most easily seen in the free 
theory where a simple computation of the scalar condensate leads to a behaviour 
(see also ref. [3]) 



oc — log 

e e 



( 



1/a 2 + m 2 



m 



log 



13 + (m/e) 2 
(m/e) 2 



) 



(47) 



e 



where we have inserted /5 = l/a 2 e 2 and m denotes the continuum fermion mass. 
Fixing m/e, as we do in this work, and approaching the continuum limit by letting 
(3 — > oo, a logarithmic increase of the scalar condensate in (3 will appear. Fig. EH 
shows clearly such an behaviour for all actions we have employed. 
It is therefore natural to use a fit ansatz of the form 



In fig. EI] we show also the fit to the data using eq. (|48|). For z = 0.2 we could 
not extract reliable values for S at f3 — 1, 2 since the determinant induced very 
large fluctuations. Since we were then left with only three data points for a three 
parameter fit, and we were also not sure whether the values for £ are possibly 
affected by finite size effects, we did not use the fit in eq.(47) for overlap fermions 
at z = 0.2. 

As can be seen, this fit function provides a nice description of the numerical 
data. In principle, for the cases of twisted mass and overlap fermions the divergent 
piece can be subtracted of from the evaluation of this term in the free theory when 
the fermion mass is matched. This is, however, very difficult for the cases of Wilson 
and hypercube fermions since there the matching of the fermion mass is not unique. 
We give in tabled the fit results using eq. (}4*S|) . 

5 Conclusions 

In this paper we have tested four different lattice fermions in their approach to the 
continuum limit in the 2-dimensional massive Schwinger model with Nf = 2 flavours 
of dynamical fermions. At fixed scaling variable z = (m^y 7 /?) 2 / 3 = 0.2,0.4,0.8 we 
have computed the pseudo scalar mass m^vT^ and the scalar condensate E for various 
values of (3 = l/e 2 a 2 . 

For all kind of fermions used, Wilson, hypercube, twisted mass and overlap 
fermions, the scaling behaviour of m^vT^ appears to be linear in a 2 . While this is 
expected for twisted mass fermions at full twist as realized here and overlap fermions, 
this result is somewhat surprising for Wilson and hypercube fermions where a linear 
dependence on the lattice spacing was expected. Of course, it might be that another 

*In 4 dimensions, there is also a logarithmic divergence. But the origin is different and can be 
removed renormalising the scalar condensate by a multiplicative rcnormalisation factor. 




(48) 
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E = A + f + Clog(^) 










(m^) 2 ' 3 = 0.2 


A 


B 


C 




Wilson 


0.118(103) 


0.215(150) 


-0.0397(447) 


0.56 


HYP 


0.0938(898) 


0.259(134) 


-0.0610(389) 


0.52 


TM 


0.0352(1189) 


0.371(179) 


-0.0767(518) 


1.7 


(m^) 2 ' 3 = 0.4 


A 


B 


C 


xVdof 


Wilson 


0.258(52) 


0.107(76) 


-0.0396(229) 


0.67 


HYP 


0.316(38) 


0.0790(543) 


-0.0357(174) 


0.48 


TM 


0.218(69) 


0.293(98) 


-0.0698(307) 


0.47 


OV 


0.249(19) 


-0.0565(271) 


-0.0200(83) 


0.13 


(m^) 2 / 3 = 0.8 


A 


B 


C 


X 2 /dof 


Wilson 


0.235(14) 


0.129(21) 


-0.145(6) 


0.62 


HYP 


0.311(14) 


0.0843(207) 


-0.162(6) 


0.19 


TM 


0.465(31) 


0.178(45) 


-0.106(14) 


2.4 


OV 


0.128(10) 


0.136(18) 


-0.160(4) 


0.022 



Table 6: Fit results for the scalar condensate using eq. (|4*%|) 

quantity can show different lattice artefacts and the scaling behaviour shows another 
dependence on the lattice spacing. 

In fig. ^2 we show our final results for m n \/j3 computed non-perturbatively and 
extrapolated to the continuum such that a direct comparison to analytical predic- 
tions can be made. Only for a value of z = 0.2 there seems to be a consistency 
with eq. (|35[) valid for strong couplings and small masses while at z — 0.8 eq. (|37[). 
obtained from a large mass expansion, seems to describe the data. In general our 
conclusion is that to the precision we could compute our results here, the analytical 
formulae do not describe the non-perturbatively obtained values of m^y^ satisfac- 
tory at all values of z. 

As a second quantity we looked at the scalar condensate. We demonstrated that 
in our 2-dimensional setup the use of a subtracted scalar condensate as derived from 
the integrated Ward identity is very useful to compute the scalar condensate since 
the so defined scalar condensate is free of divergence terms oc 1/a. Our data are also 
consistent with a logarithmic divergence in a as can be derived in the free theory. 

We also discussed some attempts to simulate the overlap operator dynamically 
by using an infrared safe kernel to construct an approximate overlap operator. This 
approximate overlap operator is then used in the simulation and physical observ- 
ables are corrected by reweighting with the determinant ratio of the exact to the 
approximate operator. We tested this idea by computing stochastically this determi- 
nant ratio. For our best candidate for an approximate overlap operator in eq. (JHBJ), 
for which we use an explicit infrared regulator in the sign function, we found that 
the eigenvalue spectra of the exact and approximate operators are very similar, see 
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fig. which shows the lemon shaped difference spectra for various accuracies of the 
approximation corresponding to different values of S in eq. (}3*3*j) . Nevertheless, even 
in this case, the fluctuations in the determinant ratio appeared to be very large even 
for a small value of the parameter 5. 

Although this led us to conclude that this stochastic way of incorporating the 
correction of the determinant ratio is not successful, we are still exploring to use 
the approximate overlap operator as the guidance Hamiltonian in the molecular 
dynamics part of the Hybrid Monte Carlo simulation while for the accept/reject 
Hamiltonians the exact overlap operator is used. We tested this setup employing 
the hypercube operator as an approximate overlap operator as the guidance Hamil- 
tonian in the molecular dynamics part. We found that even with this operator, 
which appeared to be the worst choice in the case of the stochastic estimate of the 
determinant ratio tested here, the acceptance rates look reasonable, see ref. |33j . 
We are investigating this promising result further at the moment with the other 
approximate overlap operators discussed in this paper. 
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